Transverse confinement of waves in 3D random media 
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We study the transmission of a tightly focused beam through a thick slab of 3D disordered 
medium in the Anderson localized regime. We show that the transverse profile of the transmitted 
beam exhibits clear signatures of Anderson localization and that its mean square width provides a 
direct measure of the localization length. For a short incident pulse, the width is independent of 
absorption. 



Anderson localization is a general wave phenomenon 
that results from interferences of multiply scattered 
waves in random media and leads to exponential suppres- 
sion of wave propagation beyond the localization length £ 
[ESS, Si- Recently, important advances have been made 
towards the observation and the theoretical description 
of Anderson localization of classical and, in particular, 
electromagnetic waves [IHHBIlllSOllIIll- Si S~ 
natures of localization can be found in the thickness de- 
pendence of the average transmission coefficient of a ran- 
dom sample [5 L in the shape of the coherent backscat- 
tering cone [a, [Z| , in the statistical distribution of trans- 
mitted intensity M, or in the time-of-flight profiles of 
transmitted waves [3, 13, U, 13] • In a recent experiment, 
Schwarz et al. [l]| have studied the transverse distribu- 
tion of intensity of a Gaussian light beam in a medium 
that is translationally invariant along the direction z of 
beam propagation, but disordered in the perpendicular 
xy plane. As was predicted by de Raedt et al. [3], the 
variable z plays the role of time in this case and the dis- 
tribution of intensity in the xy plane behaves very much 
as the intensity of an expanding wavepacket in a two- 
dimensional (2D) disordered medium. Because all waves 
are localized in 2D [l5j |. the expansion of the beam is 
halted after one localization length £. This phenomenon 
was called "transverse localization" 
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Even though the concept of transverse localization im- 
plies 2D disorder which is translationally invariant along 
the direction of beam propagation, one can formally at- 
tempt the same experiment in a medium with 3D dis- 
order. As has been shown very recently, the transverse, 
2D profile of intensity of a wave transmitted through a 
disordered slab shows clear signatures of genuine Ander- 
son localization in 3D [l|| • The purpose of this Letter is 
to present a detailed theoretical study of transverse con- 
finement of waves due to Anderson localization in 3D. 
This provides a guide for future experiments to exploit 
this interesting new phenomenon to access localization 
of various types of waves (light, sound, etc.). We first 
consider a continuous monochromatic beam focused to 
a point at the surface of a disordered slab and show 
that the mean square width a 2 of the spatially-resolved 
transmission coefficient T(p) scales as L£ in the local- 
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FIG. 1: A beam is focused on the surface of 3D disordered slab 
of thickness L. The bell-shaped average position-dependent 
transmission coefficient T(p) has a mean square width a 2 that 
depends on the strength of disorder in the medium. 



ized regime, when £ <C L, in contrast to the relation 
a 2 oc L 2 in the diffuse regime and even at the mobil- 
ity edge. This result suggests that a measurement of a 2 
gives direct access to the localization length £. However, 
absorption appears to produce the same effect on a 2 as 
localization, i.e. in the diffuse regime we find a 2 oc LL a , 
if the macroscopic absorption length L a is much smaller 
than L. We then replace the continuous beam by a short 
pulse. In the localized regime, the mean square width 
cr 2 (t) of the time-dependent, spatially resolved transmis- 
sion coefficient T(p, t) first grows roughly as t 1 / 2 and 
then converges to a constant value a 2 ^ oc L£, in sharp 
contrast to the relation cr 2 (t) oc Dt valid in the diffuse 
regime. Again, the localization length £ can be extracted 
directly from <J 2 (t), but contrary to the stationary case, 
<J 2 {t) does not depend on absorption. Finally, exactly at 
the mobility edge, cr 2 (t) initially grows as t 1 / 2 as in the 
localized regime, until it saturates at a 2 ^ ~ L 2 in the 
long-time limit. 

A beam focused to a point p — at the surface z = 
of a disordered slab produces a bell-shaped intensity pro- 
file T(p) at the opposite face z = L of the slab, as we 
schematically show in Fig. [TJ We analyze T(p) in the 
framework on the self-consistent (SC) theory of local- 
ization [171, extended to account for finite-size effects 
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TABLE I: Summary of analytical results for the position- 
dependent diffusion coefficient D(z, 0) and its associated 
mean square width a 2 of the position-resolved transmission 
coefficient T(p) describing transmission of a tightly focused 
monochromatic beam through a disordered slab of thickness 
L>f is the mean free path, k is the wave number, z c ~ 3£, 
£ is the localization length, and z = min(z, L~z). We assume 
z c < L for kl = 1 and I < £ < L for k£ < 1. 



0, 0, El- The intensity Green's function C(r,r',0) 
obeys a diffusion equation [— VDV]C = S(r-r') with 
boundary conditions C =F 2 o(-D / Db)9 z C — at the sur- 
faces z = and z = L of the slab and a self-consistently 
determined, position-dependent diffusion coefficient that 
obeys 



(1) 



Here k is the wavenumber, £ is the mean free path, 
Db is the Boltzmann diffusion coefficient, and zq ac- 
counts for internal reflections at the sample surface. The 
position-dependent diffusion coefficient D(r, ft) arises 
naturally from the microscopic derivation of SC theory 
in a bounded medium [l9j |. in contrast to the previ- 
ously employed concept of scale-dependent D = Z?(q, fl) 
(see, e.g., Refs. and [IH). I n the s ^ &1 ° geometry, 
D(r, Q) = D(z, f2) due to the translational invariance 
of the return probability C(r, r, f2) in the xy plane. T(p) 
is found as T{p) = -D(z = L, 0)<9 z C*({p, z = L}, {p' = 
0, z 1 = £},0). We now focus on the mean square width 
of the transmitted beam: a 2 = J p 2 T(p)d 2 p/ J T(p)d 2 p 
and consider three different regimes of wave propagation: 
weak disorder (diffusion, k£ ^> f), strong disorder (local- 
ization, k£ < 1) and the critical regime (mobility edge, 
kl=l). 

In the limit of weak disorder (k£ ^> I) the posi- 
tion dependence of D can be neglected and one can set 
D(z,0) = D B [l - [k£Y 2 ]. Straightforward solution of 
the diffusion equation then yields 



2 

°dif 



2L 2 
3 



1 -' r 



(2) 



For strong disorder (k£ < I), our self-consistent equa- 
tions do not have an exact analytical solution. We set 



D(z,0) ~ D(0, 0) exp(-2z/£) with ~z = min(z, J L- z), 
which is inspired by the observation that this D(z, 0) rep- 
resents an asymptotically exact solution in a semi-infinite 
medium for z 3> £ fl8j . Comparison with a numerical 
solution shows that this expression constitutes a rather 
good approximation everywhere, except in the central 
part of the slab, where the real value of D is larger. The 
expression for a 2 that follows from this ansatz is 
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(3) 



At the mobility edge (k£ = I) an approximate expres- 
sion for D is D(z,0) ~ D(0,0)/(1 + z/z c ). Just like in 
the case of strong disorder, this approximation for D is 
inspired by the solution for semi-infinite medium [3] and 
deviates from the numerical solution in the central part of 
the slab only; D(0, 0) and z c depend on the exact value 
of zq in the boundary conditions. For our purposes, it 
will be sufficient to set zq — and use the corresponding 
D(0, 0) = D B and z c ~ 31 This gives 
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We note that this result is not in contradiction with the 
multifractal statistics of wavefunctions at the mobility 
edge (see Ref. [22[ for a review), even though a link be- 
tween o"me an d multifractality is beyond the scope of this 
paper. 

For comparison, we summarize the analytical results 
for <j 2 following from Eqs. (HHJ) in the limit of large L in 
Table |TJ To test our analytical results, we solve the SC 
equations of localization numerically and show the result- 
ing mean square widths a 2 as functions of slab thickness 
L in Fig. O The agreement between analytical and nu- 
merical results is satisfactory, confirming the validity of 
our analytical analysis. 

An important comment is in order. It is remarkable 
that in the case of weak disorder (k£ 3> I), o 2 — 2L 2 /3 
does not depend on D and even a very small, but finite 
and spatially uniform diffusion coefficient would lead to 
the same result for a 2 . This emphasizes the importance 
of the position dependence of D(z, 0) in the localized 
regime kl < 1 and at the mobility edge kl — 1. Note 
that this is the first time that the position dependence 
of D is absolutely vital to obtain a result that is differ- 
ent from the universal diffuse outcome. Scale-dependent 
but spatially uniform diffusion coefficients D oc 1/L and 
D oc cxp(— L/£) put forward by the scaling theory [l5| 
explain the scaling of the stationary transmission coeffi- 
cient T (x 1/L 2 at the mobility edge and T oc exp(— L/£) 
in the localized regime, respectively. However, being uni- 
form in space, these expressions for D result in exactly 
the same mean square width of the transmitted beam as 
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FIG. 2: Comparison of analytical (lines) and numerical (sym- 
bols) results for the mean square width a' 2 of the spatially 
resolved stationary transmission coefficient T(p) of a slab. 



D = Db- cr 2 = 2L 2 /3. It is the position dependence of 
D that accounts for the new expressions for a 2 that we 
find at strong disorder. Incidentally, a 2 = 3L 2 /8 at the 
mobility edge differs from the result in the diffuse regime 
by just a numerical coefficient. In contrast, in the local- 
ized regime a 2 ~ 2L£ links a to the localization length 
£, which offers an elegant way of measuring the latter 
experimentally. 

Absorption was a serious obstacle for the unambiguous 
interpretation of a number of experiments on Anderson 
localization [B|, H, 0, Hj] . It is therefore important to study 
the role of absorption in the context of transverse con- 
finement of waves in 3D. A straightforward calculation in 
the regime of weak disorder (kl 1) shows that if the 
macroscopic absorption length L a = \/l/3fj, a (where \x a 
is the absorption coefficient) is much longer than the sam- 
ple thickness L, a 2 given in Table U only acquires a small 
correction — 2L 4 /45L 2 . Hence, weak absorption is not an 
obstacle for observing transverse confinement. However, 
in the opposite case of L a <C L we obtain a 2 ~ 2LL a , 
i.e. absorption plays exactly the same role as localization 
and the resulting equation for a 2 coincides with the one 
obtained in the localized regime, with £ replaced by L a . 
This indicates that a study of stationary transverse con- 
finement cannot distinguish localization from absorption 
and hence suffers from the same drawbacks as previous 
works [li, 0,0. 

A way to overcome complications due to absorption is 
suggested by recent works [1, HI U, Hj]. The idea is to 



study the dynamics of wave propagation rather than the 
stationary transport. We adopt this idea here too and 
replace the continuous incident beam in the experiment 
depicted in Fig. Q] by a short pulse. This complicates the 
analysis considerably because no simple analytic approx- 
imation exists for D(z, CI) at arbitrary CI ^ 0, contrary 
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FIG. 3: Numerically calculated mean square width a 2 (t) of 
the time-dependent transmission coefficient T(p, t) of a disor- 
dered slab in the diffuse regime (kl = 3), at the mobility edge 
(kl = 1), and in the localized regime (L/£ = 1, 2, 4). The 
dashed lines show a 2 (t) = AD B [l - (kl)~ 2 ]t for kl = 3 and 
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aZo = 2L£(1 - i/L) for L/£ = 4. The thickness of the slab is 
L = IQQl. Time is given in units of t D = (L + 2z ) 2 /-it 2 D b ■ 
The inset shows the same results in the log-log scale, the 
dashed line is u 2 (t) tx t 1 ^ 2 . 



to the stationary O = case. One way to calculate a 2 , 
which now acquires a time dependence, is to solve the SC 
equations of localization numerically. We show the re- 
sults of this calculation in Fig. [3j The difference between 
diffuse and localized regimes is manifest in Fig. [3] The 
rise of <J 2 (t) with time is unbounded for kl > 1, whereas 
cr 2 (t) saturates at a finite value cr^. for kl < 1. The latter 
result can be understood from the following approximate 
calculation which, however, turns out to be quite ade- 
quate. Instead of using the self-consistent equation for 
D(r,Cl), let us simply set D(z,Ci) — ~iCl^ 2 , which is 
the solution of SC model in the infinite medium in the 
limit CI — > 0. A straightforward calculation then yields 
0^, = 2Z£(1 - d/L) up to the first order in £/L < 1. 
This equation, shown by a dashed horizontal line in Fig. 

falls fairly close to the numerical result at L/£ = 4, 
allowing us to conjecture that it might be a good es- 
timate of a 2 ^ in the limit of L ^> £. Our calculation 
suggests that the saturation of cr 2 (t) at a constant level 
takes place not only in the localized regime, but even at 
the mobility edge. We find the asymptotic value of a(t) 
to be (Too ~ L (see Fig. [3]). It is worthwhile to note that 
our time-dependent cr 2 (t) is not sensitive to absorption 
and that all curves of Fig. [3] remain exactly the same 
in an absorbing medium. This provides a solution to the 
long-standing issue of distinguishing localization from ab- 
sorption [fl, 0, 0, [1| , as has been recently demonstrated 
in Ref. (Tf. 

Let us now analyze the time dependence of cr 2 (i). 
First, in the diffuse regime kl ^> 1 one readily obtains 
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a 2 (t) ~ ADt, where D ~ Z3_b[1 — (H) -2 ] can be assumed 
position- independent. For ki — 3 this result is shown in 
Fig. [3] by a dashed line that is indeed very close to the 
result of the numerical calculation. To study the time 
dependence of a 2 (t) at the mobility edge and in the lo- 
calized regime, we replot the curves of the main plot, 
except the one corresponding to ki — 3, in the inset of 
Fig. [3] in a log-log scale. We clearly see that the initial 
growth of u 2 (t) with t is power-law: a 2 (t) cx t a with 
ol — 0.5. It is remarkable that this power-law growth 
is observed not only in the localized regime but at the 
mobility edge ki = 1 as well. To understand the power 
exponent a ~ 0.5, let us consider a point source of waves 
located inside a random medium of size L><, very far 
from its boundaries. Neglecting the position dependence 
of D yields the return probability [13, [23| 

where q m i n ~ 1/L accounts for the finite size of the 
medium and g max ~ 1 / i ensures convergence of the inte- 
gration. Inserting this expression into Eq. ([1]) and solv- 
ing for D(fi) at ki = 1, we find D(Q) oc (-ift) 1 / 3 for 
Q » D B /L 2 and D{Q) cx (-ifl) 1 / 2 for Q < D B /L 2 . The 
typical mean square radius (r 2 ) of the intensity profile in 
3D is then (r 2 ) cx i 2/3 for t <C L 2 /D B and (r 2 ) cx t 1 / 2 for 
t 3> L 2 /D B - If we now assume that (r 2 ) ander 2 (t) = (p 2 ) 
in transmission through a slab behave in a similar way, 
we obtain a qualitative explanation for the power-law 
scaling <7 2 {t) cx t 1 / 2 that we observe in Fig. [3] at t > to- 
Note that the above simple argument does not explain 
saturation of a 2 (t) at a constant level at larger times. 

Another interesting and unexpected feature of u 2 (t) 
that we found sufficiently deep in the localized regime is 
its nonmonotonity with time: the curve corresponding to 
L/£ = 4 in Fig. [3] has a weakly pronounced but clearly 
visible maximum at t/to — 30. Even though this non- 
monotonity is predicted by the SC model (we checked 
that the maximum is present at least for L/l — 50-200 
and L/£ = 3-6), it remains to be seen whether it is an 
artifact of the model or a real physical phenomenon. 

In conclusion, we have shown that the transverse con- 
finement of waves in 3D random media can be very useful 
for demonstrating Anderson localization and for measur- 
ing the localization length in an experiment. The depen- 
dence of the mean square width a 2 of a tightly focused 
beam transmitted through a slab of random medium on 
the slab thickness L and on time (in a pulsed experiment) 
is qualitatively different in the diffuse and in the localized 
regimes of wave propagation. This qualitative difference 
is due to the position dependence of the diffusion coeffi- 
cient induced by interferences. In a pulsed experiment, 
a 2 (t) is independent of absorption in the disordered sam- 



ple and thus measuring a 2 (t) allows to avoid the risk of 
confusing localization and absorption. We expect our re- 
sults to hold independent of the microscopic structure 
of disorder. For correlated disorder, i should be replaced 
by the bare (i.e. unaffected by macroscopic interferences) 
transport mean free path. 
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